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Abstract The bonding mechanism of ethene to a nickel 
or palladium center is studied by the density matrix 
renormalization group algorithm, the complete active 
space self consistent field method, coupled cluster the¬ 
ory, and density functional theory. Specifically, we fo¬ 
cus on the interaction between the metal atom and 
bis-ethene ligands in perpendicular and parallel orien¬ 
tations. The bonding situation in these structural iso¬ 
mers is further scrutinized using energy decomposition 
analysis and quantum information theory. Our study 
highlights the fact that when two ethene ligands are ori¬ 
ented perpendicular to each other, the complex is stabi¬ 
lized by the metal-to-ligand double-back-bonding mech¬ 
anism. Moreover, we demonstrate that nickel-ethene 
complexes feature a stronger and more covalent inter¬ 
action between the ligands and the metal center than 
palladium-ethene compounds with similar coordination 
spheres. 


Keywords (i^^-transition metals • DMRG • CASSCF • 
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Fig. 1 Schematic representation of the back-donation pro¬ 
cess. M refers to either palladium or nickel. 


1 Introduction 

(i^^-transition metals like Ni, Pd, and Pt are very ver¬ 
satile metals used in batteries, alloys, and catalysts. 
In particular in organometallic chemistry, they play an 
important role in catalytic processes like coupling reac¬ 
tions [Illllinill] and cycloaddition reactions |K1IH] . Unsat¬ 
urated organic compounds like olefins can easily form 
organometallic complexes with (i^^-transition metals by 
metal-olefin bonding. It is commonly believed that the 
so-called metal-olefin bonds are formed by the process 
of back-donation. 

The Dewar-Chatt-Duncanson model miH] is widely 
used to explain the back-donation process between met¬ 
als and olefins. Metal tZ-orbitals overlap with olefin tt*- 
orbitals, allowing electron transfer from metal (i-orbitals 
to ligand 7r*-orbitals. The electron transfer process from 
bonding metal to anti-bonding ligand orbitals reduces 
the bond order of the ligand 7r-bonds. This destabilizes 
the carbon-carbon double bond and lowers the energy 
barrier to bond cleavage isiiQiiiiiiiiiiiiiiiin] . 

Recently, we elucidated the nickel-ethene reaction 
pathway and the crucial role of metal-to-ligand back- 
donation in the metal-olefin bond-formation process [16] . 
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Our study reveals the presence of a transition state 
along the nickel-ethene reaction pathway. This pecu¬ 
liar feature in the metal-olefin bond formation pro¬ 
cess motivates this study of other d^^-transition metal 
complexes containing a nickel and palladium centers 
and their reactions with small olefin ligands jOlfTT] . The 
bond-formation process of d^^-transition metals and olefins 
can be dissected using the energy decomposition analy¬ 
sis (EDA) and an orbital entanglement analysis, which 
allows us to identify the most important orbital inter¬ 
action along the reaction coordinate. 

In the EDA developed by Morokuma HUS], Ziegler 
and Rank the quantum system is divided 

into disjoint fragments according to the interaction of 
interest. In this work, the interaction energy between 
ligand(s) and the metal center, AEint, is decomposed 
into three main components, 

Z\Eint = /\Ielstat + Z\Epauli + AEq!, (1) 

where A14istat denotes the electrostatic interaction en¬ 
ergy, AEpauii is the repulsive Pauli interaction, and 
AEoi denotes the orbital interaction between the frag¬ 
ments. The EDA has proven to be a very powerful tool 
for analysing chemical bonds and orbital interactions in 
many complex chemical systems, including transition 
metal complexes [^1^1^ . 

Quantum information theory allows us to quantify 
the interaction and correlation of orbitals and orbital 
pairs The entanglement between one or¬ 

bital and the orbital bath is measured by the von Neu¬ 
mann entropy of the reduced density matrix of the or¬ 
bital of interest, here referred to as one-orbital reduced 
density matrix. The eigenvalues of the one-orbital re¬ 
duced density matrix ujo, are used to calculate the single¬ 
orbital entropy m, 

s(i)* = - 

a 

We refer the interested reader to refs for more 

details on how to calculate orbital-reduced density ma¬ 
trices. Similarly, the entanglement of two orbitals with 
the orbital bath is quantified by the two-orbital entropy 

~ ~ ; (3) 

a 

where are the eigenvalues of the two-orbital re¬ 

duced density matrix. 

The correlation between orbital pair i and j can be 
measured by the orbital pair mutual information 

mm. 

Ii\j = - ^(l)i - ^(l)i)(l - (4) 

where Sij is the Kronecker delta. 


Both s(l)i and Ii\j quantify orbital interactions and 
can be used to identify different types of electron corre¬ 
lation effects |2^[34l[35] . dissect chemical bonding j36l 
l3Qll37[l38[l39j . and locate transition state structures in 
molecular systems ioiiiiiiii]. 

In this work, we investigate the bonding situation 
in nickel-ethene and palladium-ethene compounds us¬ 
ing wavefunction approaches such as the complete ac¬ 
tive space self-consistent field approach, the density ma¬ 
trix renormalization group algorithm, and coupled clus¬ 
ter theory. In particular, we investigate the potential 
energy surfaces resulting from the interaction of the 
ethene molecule(s) approaching the palladium center in 
three structural rearrangements Pd(C 2 H 4 ), Pd(C 2 H 4 ) 2 '\ 
and Pd(C 2 H 4 ) 2 “^\ where || indicates that the ethene 
ligands are aligned in parallel, while T indicates a per¬ 
pendicular orientation of the ethene ligands (see Eig- 
ure[2|). In case of the nickel-ethene, we augment our pre¬ 
vious analysis of the Ni(C 2 H 4 ) reaction pathway with 
the symmetric bond formation process of Ni(C 2 H 4 ) 2 ^^ 
and Ni(C 2 H 4 ) 2 “^\ Eurthermore, the bonding interaction 
between the transition metals (Ni, Pd) and the ethene 
ligand(s) is analysed in terms of the energy decompo¬ 
sition analysis as implemented in ADE and an 

orbital entanglement analysis. 

This paper is organized as follows. Section [2] con¬ 
tains the computational details. In section [3l the bond 
formation process of different nickel-ethene and palladium- 
ethene complexes is dissected using the EDA and or¬ 
bital entanglement analysis. Einally, we conclude in sec¬ 
tion [H 


2 Computational Details 

2.1 Geometry Optimization 

The structures of all metal-ethene complexes are opti¬ 
mized by scanning the nickel-carbon bond from 1.75 A 
to 2.75 A and the palladium-carbon bond from 2.05 A 
to 4.05 A (constrained geometry optimization). In ad¬ 
dition, for Ni(C 2 H 4 ) 2 , a second reaction pathway was 
investigated where one ethene molecule approached the 
Ni(C 2 H 4 ) fragment. In this asymmetric bond formation 
process, all nickel-carbon distances were fixed, while all 
hydrogen atoms were allowed to freely relax. All cal¬ 
culations were performed with the Adf 2013 software 
package [i21l4;ill44j . Scalar relativistic effects were in¬ 
cluded using the ZORA Hamiltonian jl5lH6llT7] . In all 
calculations, a DZP [JS] basis set and the BP86 jl^[5Q] 
exchange-correlation functional were used. 
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Fig. 2 Equilibrium structures of nickel- and palladium-ethene complexes optimized by BP86. The || symbol indicates that 
both ethene ligands are aligned parallel to each other, while the T symbol is used to label the perpendicular arrangement of 
the ethene ligands. 


10.00 



2.5 3.0 3.5 4.0 

Pd-C distance [A] 


(a)Pd(C2H4) 




Fig. 3 Potential energy surfaces for the palladium-ethene reaction pathway in different structural rearrangements. In case of 
the two ethene molecules, the potential energy surfaces result from the symmetric dissociation of the ethene ligands from the 
metal centre. The last point of the reaction coordinate is adjusted to zero. 


2.2 CASSCF 

CASSCF [CTf5^[53] calculations for nickel-ethene were 
performed in the Dalton 2013 [5l] software package, 
while the Molpro2012 [551156] software suite was used 
for palladium-ethene. A TZP ANO-RCC basis set was 
employed in all CASSCF calculations with the following 
contraction schemes: H:(8s4p3(il/) ^ [dsApSdlf] [57] . 
C:{8s7p4:d3f2g) [4:s3p2dlf] [58] . 

Ni:(1059p8(i6/4^2h) ^ [6sbp3d2flg] [59], and 

Pd:(2l5l8pl3d6/4^2h) ^ [10s9p9(i6/4^2h] [59]. Scalar 
relativistic effects were included by the second-order 
Douglas-Kroll-Hess Hamiltonian [60l|6l]. The CASSCF 
orbitals were visualized using the Jmoll4.2.7[62] visu¬ 
alization software. 

For Ni(C 2 H 4 )^-^^ (n = 3) as well as Pd(C 2 H 4 )^"^^ 
and Pd(C 2 H 4 ) 2 '^ (n = 4), we correlated 14 electrons 
in 14 orbitals, including the nd^y^ nd^z^ ndyz, ndz 2 , 
ndx 2 _y 2 ^ and {n^l)dxy, {n^l)dxz, {n+l)dyz, {n^l)dz 2 , 

{n + I)dx 2 -y 2 orbitals from the d^^-transition metals 
and both tt and tt* orbitals from the ethene ligands. For 
Ni(C 2 H 4 ) 2 ^^, the Ni 4da,y-orbital was excluded resulting 
in CAS(14,13)SCF calculations. 

For the Pd(C 2 H 4 ) complex, we performed CAS(12,12)- 
SCF and CAS(14,14)SCF calculations. The CAS(12,12) 


active space contains the metal 4d and the ethene tt— 
and TT*— orbitals. To evaluate the contribution of the 
(7 and cr* orbitals, the CAS(12,12) active space was ex¬ 
tended to 14 electrons correlated in 14 orbitals in our 
CAS(14,14)SCF calculations. The resulting CASSCF 
orbitals along the potential energy surfaces are pre¬ 
sented in Figures S1-S18 of the Supporting Informa¬ 
tion. 

C 2 v symmetry was imposed for Pd(C 2 H 4 ), D 2 h sym- 
metry for Ni(C 2 H 4 )^"^ and Pd(C 2 H 4 )^"t and D 2 sym¬ 
metry for Ni(C2H4)2''“^ and Pd(C 2 H 4 ) 2 ''“t 


2.3 UCCSD and UCCSD(T) 

The Unrestricted Coupled Cluster Singles Doubles (UCCSD) 
and UCCSD and perturbative Triples (UCCSD(T)) [63] 
calculations were performed with the Molpro2012 [551 
[56] program. The core orbitals were kept frozen, while 
all virtual orbitals were correlated. The same basis sets, 
point group symmetries, and relativistic Hamiltonian 
were used as in our CASSCF calculations. The UCCSD(T) 
energies are collected in Tables S1-S5 of the Supporting 
Information. 
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2.4 EDA 

The energy decomposition analysis calculations were 
performed for the nickel- and palladium-ethene com¬ 
plexes at equilibrium distance using the ADF2013 [42l 
software package. Specifically, the supra-molecule 
was divided into one fragment containing the metal cen¬ 
ter and a second fragment containing the ethene ligand 
for monoligated complexes, while a third fragment com¬ 
prising the second ethene ligand was added for the bi- 
ligated metal compounds. 


2.5 DMRG 

The Budapest DMRG [64] program was used to per¬ 
form the DMRG calculations. As orbital basis, the nat¬ 
ural orbitals obtained from the largest GASSGF cal¬ 
culations as described in the previous subsection were 
used. For the biligated nickel-complexes, the active spaces 
were extended by including additional occupied and vir¬ 
tual natural orbitals. 10 additional occupied orbitals 
(2xAg, 2xB3^, 1xB2u5 l^B^g, lxBi^,lxB2g, lxB3g 
and IxAu for Ni(G 2 H 4 ) 2 '\ and 2xA, 2xBi, 3xB3 and 
3 xB 2 for Ni(G 2 H 4 ) 2 “^^) and 10 virtual orbitals for Ni(G 2 H 4 
(^2xAg, 1xB3^, 1xB2u5 l^B^g, lxBi^,lxB2g, 2xB3g 
and lx Au) and 9 virtual orbitals for Ni(G 2 H 4 ) 2 ^^ (4x A, 
IxBi, 2 xB 3 and 2 XB 2 ) were added to the active space, 
increasing it to 34 electrons correlated in 33 orbitals 
(DMRG(34,33)). The DMRG calculations for Pd com¬ 
plexes were carried out with the same active spaces as 
in GASSGF. Furthermore, we made sure that the ac¬ 
tive spaces contained similar orbitals along the reaction 
coordinate, i.e., molecular orbitals with similar atomic 
contributions. 

To enhance convergence, we optimized the orbital 
ordering [36] . The initial guess was generated using the 
dynamically extended-active-space procedure (DEAS) [31] 
In all DMRG calculations, the Davidson diagonaliza- 
tion threshold was set to 10“^ for the nickel-complexes, 
and 10“^ for the palladium compounds. The minimum 
number of block states, m, was set to 64 (in the pre¬ 
optimization), while the maximum number was set to 
1024. The convergence of DMRG with respect to m is 
summarized in Tables S1-S3 of the Supporting Infor¬ 
mation. 

The orbital entanglement and correlations diagrams 
were determined from the DMRG wavefunctions as de¬ 
scribed in ref [30] . 




)^lg. 4 Potential energy surfaces for the nickel-ethene reac¬ 
tion pathway in different structural rearrangements. The po¬ 
tential energy surfaces correspond to the symmetric dissoci¬ 
ation of the ethene ligands from the metal centre. The last 
point of the reaction coordinate is adjusted to Zero. 

3 Numerical Results 

3.1 Geometries and Potential Energy Surfaces 

The DFT-optimized geometries along the metal-ethene 
dissociation pathway are summarized in the Supporting 
Information. Figure [2] shows the equilibrium structures 
of all investigated metal-ethene compounds. For all op¬ 
timized structures, the hydrogen atoms of the ethene 

Table 1 Bonding energies for metal-ethene complexes in 
kcal/mol for CASSCF and CCSD(T). 


Molecular 

CASSCF 

UCCSD(T) 

Ni(C2H4)^"^ 

>39.0 

n/a 

Ni(C2H4)y^ 

>39.0 

n/a 

Pd(C2H4) 

6.6 

29.0 

Pd(C2H4)^"’ 

14.8 

51.4 

Pd(C2H4)y’ 

16.8 

54.5 
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Table 2 Energy decomposition analysis for nickel- and palladium-ethene complexes, re is the equilibrium distance between 
the metal center and the carbon atom of the ethene molecule(s). zAEeistat is the electrostatic interaction energy, zAEpauii is the 
repulsive Pauli interaction, and AEo\ refers to the orbital interaction between the fragments. Eint = Z\Epauii +^V"eistat +^-^oi 
is the total interaction energy. 


Molecule 

[A] 

ASpauli [S] 

AFelstat [S] 

AEoi [^1 

TTi rkcali 

Ni(C2H4) 

1.88 

385.3 

-220.1 (47%) 

-251.7 (53%) 

-86.6 

Ni(C2H4)^"’ 

2.02 

591.9 

-313.4 (45%) 

-387.2 (55%) 

-108.6 

Ni(C2H4)^-^^ 

2.02 

595.3 

-314.5 (43%) 

-411.2 (57%) 

-133.7 

Pd(C2H4) 

2.25 

125.2 

-109.4 (58%) 

-58.2 (42%) 

-42.4 

Pd(C2H4)^"^ 

2.35 

184.4 

-165.9 (66%) 

-83.4 (34%) 

-64.9 

pd(C2H4)y^ 

2.35 

184.7 

-166.5 (65%) 

-89.6 (35%) 

-77.3 


molecule are slightly bent outside the molecular (C-C- 
H) plane of the uncoordinated ethene. 

The potential energy surfaces for the palladium- 
ethene dissociation process are displayed in Figure [3l 
while Table [T] summarizes the bonding energy, z.e., the 
energy difference between the equilibrium structure and 
the dissociation limit. In general, Pd(C 2 H 4 ) 2 ^^ has the 
largest energy to association in both CASSCF and UCC- 
SD(T) calculations. The missing dynamic electron cor¬ 
relation energy in CASSCF leads to more shallow po¬ 
tential energy well depths compared to the CC results. 
In contrast to the monoligated nickel-ethene complex m, 
ah investigated quantum chemistry methods predict no 
barrier to association along the palladium-ethene reac¬ 
tion coordinate. 

Figure m shows the potential energy surfaces for the 
symmetric nickel-ethene dissociation pathway predicted 
by CASSCF, DMRG, and UCCSD. We were unable to 
converge the constrained geometry optimization for Ni- 
C distances larger than 2.8 A. Thus, only estimated po¬ 
tential well depths and bonding energies of Ni(C 2 H 4)2 
and Ni(C 2114 ) 2 ^^ are provided. Furthermore, we en¬ 
countered convergence difficulties in our UCCSD cal¬ 
culations for stretched distances of Ni(C 2 H 4 ) 2 . There¬ 
fore, it remains unclear if the symmetric dissociation of 
Ni(C 2 H 4)2 features a transition state as found for the 
monoligated nickel-ethene complex. 

As shown in Table [U the lower bound for the bond¬ 
ing energy is considerably larger in nickel-ethene than 
in palladium-ethene complexes suggesting a stronger 
bonding interaction between the nickel center and the 
ethene ligands in terms of 7r-donation and metal-to- 
ligand back-bonding. 

In general, CC calculations predict shorter metal- 
ethene bond lengths than found in CASSCF, which can 
be attributed to the missing dynamic electron correla¬ 
tion effects in the latter. Specfficahy, the CASSCF Pd- 
C equilibrium distance in Pd(C 2 H 4 ) is approximately 
2.25 A, which reduces to 2.10 A in UCCSD(T). For 
Pd(C 2 H 4 ) 2 the equilibrium bond length decreases from 


2.35 A in CASSCF to 2.25 A in UCCSD(T), while 
the equilibrium bond lengths are slightly shorter for 
Pd(C 2 H 4 )^"^t 2.35 A in CASSCF and 2.15 A in UCCSD(T). 


3.2 Elucidating the Metal-to-Ligand Back-Donation 

3.2.1 Energy Decomposition Analysis 

The EDA results at the equilibrium distance for ah 
metal-ethene complexes are summarized in Table [2j The 
total interaction energy (see Table [2]) is dehned as the 
sum of AEpauii, AWistat, and AE^i and quantihes the 
interaction between the fragments. 

Ah investigated nickel-ethene complexes have a con¬ 
siderably larger total interaction energy (in absolute 
value) than the corresponding palladium-ethene com¬ 
pounds with similar coordination sphere. Comparing 
AWistat with AEoi, we observe that AVoi constitutes 
the dominant contribution in nickel-ethene complexes, 
while Z\Wistat dominates in palladium-ethene complexes. 
Since Z\Wistat corresponds to the classic electrostatic 
interaction between fragments and AE^i represents the 
interaction between orbitals on one fragment with the 
orbitals on the other fragment, a larger contribution of 
AWistat indicates more ionic interactions between the 
fragments, while a larger contribution of AE^ suggests 
a stronger covalent nature of the interaction between 
the fragments. The different ratios between Z\Wistat and 
AEoi suggest that the nickel-ligand bond is more cova¬ 
lent, while the palladium-ligand bond is more ionic. 

Table 3 Relation between the strength of orbital entangle¬ 
ment and correlation and electron correlation effects. 


Correlation effects 

s(l)i 


Nondynamic 

>0.5 

w 10-1 

Static 

0.5-0.1 

Ri 10-2 

Dynamic 

<0.1 

w 10-3 
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Mutual information 


single-orbital entropy 
(a) =2.25A 
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Fig. 5 Orbital-pair mutual information and single-orbital entropy for Pd(C 2 H 4 ) determined from DMRG(14,14) calculations. 


3.2.2 Orbital Entanglement 

An orbital entanglement analysis uses the single-orbital 
entropy 8(1)^ to measure orbital entanglement and the 
orbital-pair mutual information I^j to quantify the cor¬ 
relation between orbital pairs. Both s{l)i and Ii\j are 
represented using diagrams. Specifically, the strength 
of the orbital-pair mutual information is colour-coded. 
Strongly correlated orbital pairs are connected by blue 


lines ^ 10“^), moderately correlated orbitals by 
red lines {Ii\j ~ 10“^), while weakly correlated orbitals 
are indicated by green lines (Ab ~ 10“^), etc. As pre¬ 
sented in ref [29], the strength of orbital entanglement 
and correlation can be associated with electron correla¬ 
tion effects (see Table [3|). Since we are interested 

in bond formation processes, our analysis will focus on 
orbitals and orbital pairs with moderately to strongly 
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Fig. 6 Orbital-pair mutual information and single-orbital entropy for Pd(C 2 H 4)2 determined from DMRG(14,14) calculations. 


entangled orbitals, z.e., orbitals with ^(l)^ > 0.1 and 
Ii\j > 10 - 2 . 

Orbital Entanglement and Correlation in Palladium- 
Ethene Figures [5ll7| show the mutual information and 
single-orbital entropy for palladium-ethene complexes 
at different points of the reaction coordinate (see Fig¬ 
ure El along with the strongly entangled molecular or¬ 
bitals. 


For the monoligated Pd(C 2 H 4 ) in the dissociation 
limit (Figure [5jc)), molecular orbitals centered on the 
metal atom (Pd Ad- and 5d-orbitals) and on the ethene 
fragment (tt- and 7r*-orbitals) are correlated. No signif¬ 
icant orbital correlations can be observed between or¬ 
bitals centered on different fragments. The most strongly 
correlated orbitals are the ethene tt- and tt* -orbitals 
(nos. 3 and 9 in Figure [5jc)). When the ethene molecule 
approaches the metal center (see Figure [5fb)), the Pd 
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Fig. T Orbital-pair mutual information and single-orbital entropy for Pd(C 2 H 4 )^ determined from DMRG(14,14) calculations. 


dyz-oihitsl (no. 8) and the ethene 7r*-orbital (no. 9) be¬ 
come weakly correlated. This orbital correlation corre¬ 
sponds to the metal-to-ligand 

back-donation process. However, the most strongly cor¬ 
related orbitals remain centered on the ethene ligand 
(tt-tt*) and the metal atom (4d-5(i), respectively. Around 
the equilibrium structure, the molecular orbitals in¬ 
volved in metal-to-ligand back-bonding (nos. 8 and 9 
in Figure [5fa)) become moderately correlated, while 


molecular orbitals involved in 7r-donation from the ethene 
TT-orbitals to the Pd 4(i;22-orbital (nos. 2 and 3 in Fig¬ 
ure [SJ^a)) are only weakly correlated. The dominant or¬ 
bital correlations remain between the ethene tt- and tt*- 
orbitals and between Pd M- and StZ-orbitals. 

A similar trend in the orbital correlation and entan¬ 
glement diagrams can be observed for Pd(C 2 H 4 ) 2 (see 
Figure [6]). In the dissociation limit (Figure (U^c)), or¬ 
bital correlations remain distributed among the ethene 
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Fig. 8 Orbital-pair mutual information and single-orbital entropy for Ni(C 2 H 4 )| determined from DMRG(34,33) calculations. 


TT- and 7r*-orbitals (nos. 3, 6, 9, and 13) as well as Pd 
4d- and 5d-orbitals (nos. 1, 2, 4, 5, 7, 8, 10, 11, 12, 
14). When both ethene ligands approach the Pd center 
(see Figure [6jb)), the Pd d^;s-orbital (no. 12) and the 
ligand (tt^ +7r|)-orbital (no. 13) are weakly correlated. 
These orbitals are involved in the metal-to-ligand back- 
bonding process. However, the dominant orbital corre¬ 
lations remain centered on the ligand orbitals and on 
the metal orbitals, respectively. Around the equilibrium 


structure, the changes in the correlation and entangle¬ 
ment patterns are more profound than for the mono- 
ligated palladium-ethene complex. While the Pd dy^- 
orbital (no. 12) and the ligand (tti+ 7r2)-orbital (no. 13) 
are moderately correlated, the correlation strength be¬ 
tween the ligand (tt^ + 7r|)-orbital and the remaining 
bonding and antibonding combinations of the tt- and 
7r*-orbitals decreases. Similarly, the Pd (i^ 2 -orbital and 
the ligand (tti +7r2)-orbital are weakly entangled,! sug- 
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Fig. 9 Orbital-pair mutual information and single-orbital entropy for Ni(C 2 H 4 )^ determined from DMRG(34,33) calculations. 


gesting a negligible contribution of 7r-donation in the 
bond-formation process of Pd(C 2 H 4 ) 2 '^ 

For Pd(C 2 H 4 ) 2 ^\ the Pd d-orbitals and the ligand 
TT- and TT*-orbitals remain uncorrelated in the dissocia¬ 
tion limit and for stretched palladium-ethene distances 
(see Figure [ll^b) and (c)). In contrast to Pd(C 2 H 4 ) 
and Pd(C 2 H 4 ) 2 '\ the correlation and entanglement di¬ 
agrams drastically change around the equilibrium dis¬ 
tance (see Figure [Tfa)). At this point, the Pd Myz- 
orbital (no. 9) and the ligand (Trjf, Tr^)-orbital (no. 10) 
as well as the Pd 4(ia,;s-orbital (no. 6 ) and the ligand 
(ttJ,TT^)- orbital (no. 7) are moderately correlated. Fur¬ 
thermore, the correlation between the ligand tt and tt*- 
orbitals (nos. 3, 7, 10, and 13) reduces and all lig¬ 


and orbitals are only moderately correlated compared 
to the ligand orbitals in Pd(C 2 H 4 ) and Pd(C 2 H 4 ) 2 '^ 
These dominant correlations between metal and lig¬ 
and orbitals suggest that the electronic structure of 
Pd(C 2 H 4 ) 2 ^^ features two metal-to-ligand back-bonding 
interactions. This double-back-bonding mechanism may 
lead to an additional stabilization of the Pd(C 2 H 4 ) 2 ^^ 
isomer compared to the Pd(C 2 H 4 )P^ complex and elu¬ 
cidates the larger orbital interaction energy AE^i of 
Pd(C 2 H 4 )y^ in the EDA. 

Orbital Entanglement and Correlation in Nickel-Ethene 
The entanglement and correlation diagrams for Ni( 02114)2 
at different points along the reaction coordinate are 
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shown in Figure [8l In the dissociation limit, the leading 
orbital correlation is found between the Ni and 

45 -orbitals (nos. 4 and 6 in Figure [H^c)). In contrast 
to Pd(C 2 H 4 ) 2 '\ the metal Myz and ligand (Trjf + ttI)- 
orbitals are already moderately correlated. Their cor¬ 
relation further increases when the ethene ligands ap¬ 
proach the metal center (see Figure [ 8 fb)). Around the 
equilibrium structure, the strong correlation between 
the Ni and 45 -orbitals diminishes and the bond¬ 

ing and antibonding combination of the Ni 3dz2-OThit8i\ 
and the ligand (tt^ + 7 r|)-orbital (nos. 3 and 6 in Fig¬ 
ure [HKa)) are strongly correlated. The latter orbital cor¬ 
relation corresponds to the metal-to-ligand 7r-donation 
mechanism. In contrast to the monoligated Ni(C 2 H 4 ) 
complex (see ref [16] for details), 7 r-donation does not 
commence until close to the equilibrium geometry and 
the corresponding orbital correlations are comparable 
to those between the Ni 3dyz and ligand (ttJ + ttJ)- 
orbitals. 

The reaction pathway of Ni( 02114 ) 2 ^^ features a no¬ 
tably different evolution of orbital correlation and en¬ 
tanglement compared to its structural isomer Ni(C 2 H 4)2 
In the vicinity of dissociation (Figure[9fc)), the Ni 8 ^ 3,2 _^ 2 - 
and 3 d; 22 -orbitals (nos. 4, 6 , and 7) as well as the lig¬ 
and (tti, 7 r 2 )*-orbital (no. 5) are strongly correlated with 
each other, while the orbitals involved in metal-to-ligand 
back-bonding (nos. 21 and 22 as well as nos. 29 and 
30) are moderately correlated. When the ethene ligands 
approach the metal center (Figure [9jb)), metal and 
ligand orbitals that participate in 7 r-donation (nos. 5 
and 6 ) and metal-to-ligand back-bonding (nos. 21 and 
22 as well as nos. 29 and 30) are strongly correlated. 
Close to the equilibrium geometry (Figure [9fa)), we 
observe a transition of orbital entanglement and corre¬ 
lation patters. Specifically, the correlation between the 
Ni 3d- and 4d-orbitals (nos. 15 and 16 as well as nos. 4 
and 7) decreases when approaching the equilibrium ge¬ 
ometry. Similar to Pd(C 2 H 4 )^, the orbital correlation 
and entanglement analysis suggest that Ni(C 2114 ) 2 ^^ 
features two metal-to-ligand back-bonding interactions 
which may stabilize the Ni(C 2114 ) 2 ^^ complex compared 
to the Ni(C 2 H 4 ) 2 '^ isomer. 

Comparison of the Bonding Mechanism in Ni- and Pd- 
Olefines Finally, we will compare the bond formation 
process and the bonding interactions in nickel-ethene 
and palladium-ethene complexes. Both the EDA and 
orbital entanglement analysis highlight the different na¬ 
ture of the metal-ethene bond and of the bond-formation 
mechanism. In general, the bonding interaction in the 
nickel-ethene complexes is stronger than in the corre¬ 
sponding palladium-ethene compounds. Furthermore, 
the degree of covalency of the metal-olefin bond is higher 


for nickel-ethene than for palladium-ethene ( cf. the large 
values of AE^i for Ni(C 2114 ) 3 , compared to Pd( 02114 ) 3 , 

{x = 1,2). The different bonding nature and bond- 
formation processes is supported by our orbital entan¬ 
glement analysis. While the metal-to-ligand back-bonding 
mechanism plays an important role in the bond-formation 
process in nickel-olefin compounds, and which estab¬ 
lishes for stretched Ni-C 2 H 4 distances in the vicinity of 
dissociation, the metal-to-ligand back-bonding in Pd(C 2 H 4 ) 3 , 
becomes important close to the equilibrium geometry. 
Similarly, the role of 7 r-donation considerably differs in 
nickel- and palladium-ethene. Specifically, our entan¬ 
glement analysis predicts that the 7r-donation mecha¬ 
nism is insignificant in Pd(C 2114 ) 3 , complexes, while it 
forms an essential part in the bond-formation process in 
Ni(C 2 H 4 ) 3 , compounds where the correlation between 
the ligand 7r-orbitals and the metal (i; 22 -orbital increases 
when the ethene ligands approach the nickel center. 

4 Conclusions 

In this work, we investigated the interactions between 
ethene ligands and the nickel and palladium center along 
the metal-ethene reaction coordinate and for perpen¬ 
dicular and parallel orientations of the ethene ligands. 

While both nickel and palladium are d^^-transition met¬ 
als, they exhibit a considerably distinct bonding mech¬ 
anism and interactions with ethene ligands. Specifically, 
nickel-carbon bonds are shorter and stronger than palladium- 
carbon bonds for both the parallel and perpendicu¬ 
lar orientation. Moreover, the bond between nickel and 
ethene has predominantly covalent character, while the 
palladium-ethene bond has mainly ionic character. 

Both d^^-transition metals create more stable com¬ 
plexes with ethene in perpendicular orientation, where 
two metal-to-ligand back-bonding mechanisms can be 
observed. The double-back-bonding allows for stronger 
orbital interactions. Moreover, our entanglement anal¬ 
ysis indicates that molecular orbitals involved in tt- 
donation from the ethene 7r-orbitals to the metal d ^2 - 
orbital are considerably more correlated in nickel-ethene 
than palladium-ethene complexes. Thus, while 7 r-donation 
plays an important role in the bond-formation process 
of nickel-ethene, the palladium-ethene bond does not 
feature strong 7r-donation. 

This work demonstrates that concepts from quan¬ 
tum information theory constitute a useful and com¬ 
plementary tool to well-established methods like energy 
decomposition analysis in dissecting chemical reactions. 
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